The Hubble Deep Field in the Far Ultraviolet 



Henry C. Ferguson 

Space Telescope Science Institute 
3700 San Martin Drive 
Baltimore, MD 21218 USA 



Abstract. Results from a recent HST survey of field galaxies at wavelengths 
1600 A and 2400 A are be presented. The data are used to constrain the fraction 
of Lyman-continuum radiation that escapes from galaxies at redshifts z fa 1. The 
combined UV-IR photometry for HDF galaxies is also used to investigate whether 
low-mass starburst galaxies dominate the field-galaxy population at redshift z ~ 1. 
The relative lack of objects with the colors of faded bursts suggests that star- 
formation is largly quiescent rather than bursty or episodic. 



When the ultimate self-consistent model of cosmology and galaxy evolu- 
tion comes along it will be able to reproduce the full distribution of galaxy 
shapes, sizes, masses, and spectral types, at all redshifts and in all envi- 
ronments. Until that happens, we as observers are forced to try to piece 
together portions of the puzzle, gathering data as new facilities and tech- 
niques become available, and testing and refining parts of the theory. One of 
the new facilities that has relatively recently become available is the Space 
Telescope Imaging Spectrograph (STIS), which was installed on the Hubble 
Space Telescope in 1997. Prior to STIS, the deepest observations of galaxies 
at far-UV (A <; 2500 A) wavelengths reached to m ~ 20.7. [] Far-UV obser- 
vations provide a sensitive probe of the population of young massive stars 
in galaxies, and spectral breaks at Lyman a and the Lyman break (912 A) 
are among the strongest features that can be used for estimating photomet- 
ric redshifts. The STIS Multi-Anode-Microchannel Array (MAMA) detectors 
provide sensitive imaging over a 25" x 25" in two bands centered at approx- 
imately 1600 A (FUV) and 2400 A (NUV). Rest-frame emission shortward 
of 2000 A can be measured for the first time in galaxies in the crucial red- 
shift range 0.5 < z < 2.5, probing the period where most of the stars in the 
present-day universe were formed. 

Amongst the aims of our survey were to (1) measure galaxy counts in the 
far-UV, (2) compare rest-frame UV morphologies to optical morphologies for 
galaxies at 0.5 <J z <; 2, (3) measure the far-UV background and constrain the 
contribution from faint galaxies, (4) constrain the contribution of galaxies to 
the metagalactic ionizing background radiation, and (5) constrain the star- 
formation histories of dwarf galaxies at redshifts z < 1. The observations 



1 Magnitudes throughout this contribution are on the AB system (Oke 1974), where 
AB = -2.5 log /„(njy) +31.4 
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cover an area of one square arcminute, with three sigma detection limits 
of AB — 29.1 in both the NUV and FUV bands. Galaxy counts from our 
survey have been presented by Gardner et al. JTcfl and constraints on the 
FUV background by Brown et al. ||. In this contribution, we concentrate 
on two topics: setting constraints on the ionizing radiation escaping from 
galaxies and setting constraints on star-formation histories of low-luminosity 
galaxies in the HDF. 

1 The Ionizing Radiation Escaping from Galaxies 

Studies of quasar absorption lines indicate that the intergalactic medium 
is highly ionized out to redshifts of order z = 5. Instead of appearing as a 
smooth Gunn-Peterson (1965) absorption trough at redshifts less than 1216 A 
in the rest frame QSO of the QSOs, the absorption due to neutral hydrogen in 
the IGM appears at discrete wavelengths (the Lyman-alpha forest). The total 
baryon content in neutral gas is much less than that expected from primordial 
nucleosynthesis, and both the measurements of the line strengths of different 
ionization states heavier elements, and the measurement of the proximity 
effect suggest that the gas is highly ionized by a metagalactic radiation field 
J u ~ 10~ 21 ergs^ 1 cm -2 Hz -1 sr _1 E3. It is likely that the ionization of 
the IGM plays an important role in regulating the cooling and collapse of gas 
to form galaxies JT]j3~S|] - 

The IGM at high redshifts is most likely ionized by young galaxies and 
QSOs. However, the relative contribution of these two types of objects to 
the high-rcdshift ionizing background is not well constrained. The declining 
space density of QSOs beyond z — 3 suggests that even if they produce 
most of the ionizing photons at z ~ 2, their contribution is likely to decline 
dramatically at high redshifts. An extrapolation of the measured luminosity 
functions and redshift evolution suggests QSOs and other AGN are likely 
to provide only about 25% of the ionizing photons at 1 Ryd at z=5. In 
contrast, the Lyman continuum flux estimated to arise from the integrated 
LF of galaxies at z ~ 3 appears sufficient to keep the IGM ionized, provided 
that a fraction / csc > 50% of the ionizing photons escape 

The most direct constraints on / esc at redshift z — come from observa- 
tions of four starburst galaxies with the Hopkins Ultraviolet Telescope (HUT) 
[ pT| . More recently Steidel et al. |l7j have reported a detection of Lyman con- 
tinuum radiation from a composite spectrum of 29 Lyman-break galaxies at 
redshifts z = 3.40 ± 0.09. The measured ratio of 900 A to 1500 A emission 
-Fgyxg is remarkably high, indeed higher than expected from some popula- 
tion synthesis models even if 100% of the of the intrinsic Lyman continuum 
radiation escapes. However, the measurement of this break in the composite 
spectrum is extraordinarily difficult, and a substantial correction for IGM at- 
tenuation is required to derive the intrinsic flux decrement of the composite 
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galaxy spectrum. Even if the detection itself is real, there are concerns that 
the sample is biased toward the bluest of the Lyman-break galaxies. 

There are seven galaxies in our STIS HDF field with spectroscopic red- 
shifts 1 < z < 1.3. At z = 1.1 the center of the FUV bandpass is at 
Arest ~ 760 A, and the integrated throughput from wavelengths longer than 
Arest = 912 A is less than 0.2% of the total throughput. To maximize the 
signal-to-noise ratio, UV fluxes and upper limits for the seven galaxies have 
been estimated using a \ 2 fitting technique (TFIT) Jl5| on the STIS images 
after convolution to the resolution of the F450W HDF images. 
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Fig. 1. Photometry for an HDF galaxy at z = 1.01. The points with errorbars show 
the measured photometry. The galaxy is undetected in the shortest wavelength 
point; the la uncertainty is shown as a solid error bar, and the 2a uncertainty as 
a dotted error bar. The curves show a fiducial model SED (not fit to the data) for 
an unobscured galaxy of age 10 9 yr undergoing constant star formation with solar 
metallicity, and a Salpeter IMF. The top curve is the emergent spectrum from the 
galaxy; the bottom curve is attenuated by the IGM. 



Figure 1 shows the photometry for one of the galaxies, superimposed on 
a typical star-forming galaxy spectral energy distribution with and without 
Madau IGM attenuation. Only one of the seven galaxies in our sample 
is detected, and that one at only the 2.5a level. 

For a galaxy at z — 1, the half- maximum points of the FUV bandpass are 
at rest-frame wavelengths 730 and 830 A, hence a small correction is necessary 
to estimate F(900). We have so far done this only in a crude way, accounting 
for the typical IGM attenuation but assuming a flat spectrum in /„ rather 
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than a model stellar population. The resulting values of -F9/15, displayed in 
Figure 2, should thus be regarded as preliminary, but are unlikely to change 
by more than a few percent. 
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Fig. 2. Ratio Fg/ 15 vs. UV spectral slope /3 (fx oc A' 3 ), for the HDF galaxies (solid 
triangles) compared to the results from Leitherer et al. Jll| (open triangles) and 
Steidel et al. 1l7|] (solid circle). The value of (3 has not been measured for the Steidel 
et al. composite galaxy, and it has been placed arbitrarily at /3 = 1.2. For the HDF 
sample /? is estimated from the photometry at A rcst < 2500 A, while for the Leitherer 
et al. sample it is measured from the HUT spectra in absorption-line free regions 
with 1220 < A < 1800. Nearly all of the upper limits for galaxies observed with 
HST or HUT fall below the Steidel et al. Q measurement. 



For unreddened star-forming objects, population synthesis models predict 
F9/15 ~ 0-2, which is roughly the minimum level of Lyman continuum jump 
consistent with the Leitherer et al. spectra. Leitherer et al. used Hcv 
measurements of the galaxies in their sample to set a lower limit on the total 
number of Lyman-continuum photons that must be present. Based on that 
analysis they derived upper limts on / esc ranging from 1% to 15%. Lacking 
Ha fluxes for the HDF galaxies, we must adopt a model spectrum to infer / esc 
from -F9/15. For this preliminary analysis we adopt the model shown in Fig. 
1; in the final analysis this will be refined by fitting models to each individual 
galaxy. For this model, if we assume that all the UV radiation produced by 
stars at rest-frame 1500 A escapes from the program galaxies, then our upper 
limits translate into / esc < 0.2 to 0.9 (depending on the galaxy). If we instead 
make the more reasonable assumption that there is some dust attenuation 
of the stellar radiation at both 900 A and 1500 A, then adopting a factor of 
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4.7 attenuation at 1500 A |18| and a Calzetti [|j attenuation curve leads to 
upper limits / csc < 0.05 to 0.2, significantly below the claimed detection ]T7| , 
but compatible with the expectations from recent radiative transfer models 
[Q. In general, accounting for dust and gas attenuation, / csc < 20% for all 
galaxies with sufficiently deep measurements from both HUT and HST. 



2 Star-formation in galaxies at z ~ 1 

It is likely that the ionizing background inhibits the cooling of gas in low- 
density regions, and in particular in the halos of low-mass galaxies [p|Jl] Jl9| , [l4]| . 
This has led to suggestions that a large fraction of the faint-blue galaxies are 
dwarf galaxies undergoing starbursts at moderate redshifts. While the most 
extreme of these models are ruled out Q , it is still possible that dwarf galaxies 
undergoing some form of episodic star formation are a major contributor to 
faint-galaxy counts . 

From the combined multi-wavelength HDF observations, it is possible to 
investigate whether starbursts or episodic star formation are ubiquitous or 
rare in the faint-galaxy population. The argument goes as follows: if dwarf 
galaxies are forming stars in bursts with star-formation timescales t$f and 
with intervals between episodes Z\ibmst >> Tsf, then provided this is roughly 
a steady-state process there ought to be 10 times as many galaxies with ages 
t ~ IOtsf as there are galaxies with age tsf- And there ought to be 100 times 
as many galaxies with ages t ~ IOOtsf provided Ziiburst is greater than t. We 
can thus examine the "fading sequences" in color-magnitude or color-color 
space and see whether the galaxies are piling up in location expected if they 
are bursting and then fading. 

Figure 3 shows a color-absolute-magnitude diagram of HDF galaxies with 
photometric redshifts 0.75 < z < 1.25, spanning a range of lookback times 
from 7 to 9.2 Gyr for a cosmology with h, fltot, ^m, = 0.65, 1, 0.3, 0.7. The 
two panels of the figure show the regions that should be occupied by galaxies 
with stellar masses of 10 9 M Q for with Calzetti Q extinctions 0.02 < Eb~v < 
0.4 and ages in three ranges of 0.5 dex. The left panel shows r = 10 7 yr burst 
models and the right panel shows constant star-formation models. Now, given 
that the number of years spanned by the three regions increases by a factor 
of 10 as they move to older ages, any steady state process should have the 
galaxies bunching up at the older ages. The distribution of colors for the HDF 
galaxies is clearly inconsistent with most of them undergoing short bursts. 
Under these model assumptions the HDF galaxy colors cluster around ages 
of 10 8 yr, but there is no increase of galaxy numbers toward older ages. In 
contrast, the loci of HDF points are reasonably consistent with that expected 
from roughly constant star- formation rates. From this type of analysis, a 
fairly robust conclusion is that star-formation in faint-blue galaxies at z ~ 1 
is a relatively quiescent process. Dwarf galaxies undergoing single bursts or 
episodic star formation with a duty cycle At ;> 10 9 yr are not the dominant 
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Fig. 3. Rest-frame M(1500A)— Iab color-magnitude diagram for HDF galaxies 
with photometric redshifts 0.75 < z < 1.25. Galaxies are from an H-band selected 
catalog covering the entire 5 square arcminutes of the HDF-N and rest-frame 
magnitudes have been interpolated from the multi-band photometry. Photometric 
redshifts were determined via template fitting. Superimposed on the data are shown 
the areas that should be occupied by galaxies in three ranges of 0.5 dex in age, 
with Calzetti § extinctions 0.02 < E(B-V) < 0.4. The left panel shows starburst 
models. The stellar mass in the models is 10 9 Mq, the metallicity is 1/5 solar, and 
the star-formation e-folding timescale is r = 10 7 yr, which is roughly the timescale 
expected for supernova winds to drive the gas from a dwarf galaxy. The right panel 
shows models with a constant star-formation rate of lMQyr -1 . Thus at an age 
of 10 9 yr the models in both panels have the same mass. The basic point is that 
distribution of the HDF galaxies in this color-magnitude plane makes sense if their 
star formation is relatively constant, while it is hard to understand if their star 
formation is bursty or episodic. 



population. Alternatively, the bursting galaxies must fade more rapidly than 
expected from simple population-synthesis models, either because the IMF is 
top heavy, or because the stellar population expands when the gas is driven 
out, and the galaxies then drop below the HDF detection limit in surface 
brightness. 
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